Deviation from the Cosmological Constant or Systematic Errors? 



o 
o 

(N 

5-1 

< 

o 

u 

43 : 

Oh- 
i 

o 

(N 
> 

m 
oo 
(N 

O 

&\ 
O 



13 



Shi Qi, 1 ' 2 ' 3 ^ Tan Lu, 1 - 2 ^ and Fa- Yin Wang 4 '@ 

1 Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, China 
Joint Center for Particle, Nuclear Physics and Cosmology, 
Nanjing University - Purple Mountain Observatory, Nanjing 210093, China 
3 Kavli Institute for Theoretical Physics China, CAS, Beijing 100190, China 
■» Department of Astronomy, Nanjing University, Nanjing 210093, China 
(Dated: April 29, 2009) 

Motivated by the fact that both SNe la and GRBs seem to prefer a dark energy EOS greater 
than —1 at redshifts z > 0.5, we perform a careful investigation on this situation. We find that the 
deviation of dark energy from the cosmological constant at redshifts z > 0.5 is large enough that 
we should pay close attention to it with future observational data. Such a deviation may arise from 
some biasing systematic errors in the handling of SNe la and/or GRBs or more interestingly from 
the nature of the dark energy itself. 
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It has been ten years since the discovery of the cosmic 
acceleration [l|, Q > which is attributed to the mysterious 
component — dark energy. In addition to the cosmo- 
logical constant, a lot of dark energy models have been 
proposed to explain the cosmic acceleration (see for ex- 
ample 3]). Though the standard ACDM model fits the 
observational data well, there are also a variety of other 
dark energy models could not be ruled out due to the pre- 
cision of current data. It is therefore still a crucial issue 
that whether the dark energy is simply the cosmological 
constant or not. 

Among all kinds of observational sources, type la su- 
pernova (SN la), which has been widely used as standard 
candles, is one of the most important classes of data 
that could impose significant constraints on the nature 
of dark energy. One of the important reasons for this is 
that SNe la provide data points along the redshifts and 
could therefore recover the nature of dark energy at dif- 
ferent redshifts. However, due to the limitation of the 
redshifts of SNe la, it is difficult to study the nature of 
dark energy beyond redshift of 1.7 with SNe la. While 
gamma-ray bursts (GRBs), as the most luminous astro- 
physical events observed today, extend the redshift to 
z > 6. After calibrated, they could be used as comple- 
mentary sources to SNe la at high redshifts in cosmology 
studies (see e.g. [3]), which has recently attracted much 
attention. At present, GRBs are still not as ideal stan- 
dard candles as SNe la. The scatters of known luminosity 
relations of GRBs are still very large and they have cir- 
cularity problem due to the lack of low redshift samples. 
In spite of this, works by many authors have put it for- 
ward in their cosmic applications. For example, recent 
advances include that new luminosity relations are intro- 
duced [B|, and that model-independent calibrations are 
proposed 0, H, 0, [13, EH- Among works using combined 
data of SNe la and GRBs, it is notable in [8( that the im- 
provements on the constraints made by including GRBs 
show that the dark energy equation of state (EOS) is 



slightly shifted towards w > — 1 at redshifts z > 0.5 [271 ]. 
Since the results there are still totally consistent with 
the cosmological constant at 2a confidence level and the 
inclusion of GRBs is much preliminary (Only systematic 
errors of luminosity relations are included for simplicity) , 
we cannot draw any concrete conclusion only from that. 
However, it is interesting that, when we have more SN la 
samples, SNe la themselves also show the same trend, as 
was shown by Figure 17 in 12]. (To be fair, There are 
also other analyses with earlier SN la sets as well as with 
the Union set [l2j show the sign of the possible increase 
in the dark energy EOS, see for example [3, [3, 13. 
The advantages of the results in [f| and Figure 17 in \l2\ 
are that the redshift binned parameterization is used, 
which assumes less about the nature of the dark energy 
compared to simple parameterizations especially at high 
redshifts. See, for example, in [3 ]. where the redshift 
binned parameterization are adopted with earlier SN la 
sets, the constraints on the dark energy EOS at redshifts 
z > 0.5 are still too weak. See also Figure 11 in 17j for 
an illustration of priors imposed on the dark energy by a 
simple parameterization itself.) Because GRBs and SNe 
la are independent sources, the fact that they both seem 
to prefer a dark energy EOS greater than —1 at redshifts 
z £ 0.5 may be worth our more attention. Motivated by 
this, we perform a careful investigation on this situation 
in this paper. 

In [H, five luminosity relations are used for GRBs, i.e. 
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To avoid circularity problem, calibration parameters and 
cosmological parameters are fitted simultaneously, and in 
the calculation of x 2 , only the systematic errors of the lu- 
minosity relations are taken into account for simplicity in 
the preliminary study of the evolution of the dark energy 
EOS including GRBs (see [|| for details). In this paper, 
we include also the measurement errors of GRBs for a 
careful investigation, i.e., of ot = cr^ ca + of ys , where the 
systematic errors, <7 sys s, are derived by requiring the re- 
duced x 2 equal to 1 and the measurement errors, <r mea s, 



are given by a. 
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-y + b 2 a1 for the fitting to 

y = a + bx, 



(6) 



where x and y denote the logarithm of the luminosity 
indicators and the logarithm of luminosity or energy of 
GRBs (see Eq. ((T][5|)). For asymmetric measurement er- 
rors, the errors of the side near the line being fitted to 
are used fu\ . 

In addition to GRBs, we have used Union compila- 
tion of SNe la from BAO measurement from [ij} 
and Q m h = 0.213 ± 0.023 from We assumed the 

prior Q k = -0.014 ± 0.017 [2(| for the cosmic curva- 
ture. We adopted the redshift binned parameterization 
for the dark energy EOS, as proposed in [2l|, to estimate 
possible evolution of the dark energy. In this parame- 
terization, the redshifts are divided into several bins and 
the dark energy EOS is taken to be constant in each 
redshift bin but can vary from bin to bin. And a set 
of decorrelated EOS parameters are introduced subse- 
quently by an appropriate transformation. The evolution 
of the dark energy with respect to the redshift could be 
estimated from these decorrelated EOS parameters. In 
this paper, we divided redshifts at points z = 0.2, 0.5, 1 
and Markov chain Monte Carlo techniques are used with 
<^(10 6 ) samples generated for each result. Since current 
observational data have only very weak constraints on 
the nature of dark energy at redshifts z > 1, we focus 
our analyses on the first three redshift bins. 

Figure[T]shows the results derived from data set includ- 
ing GRBs. As stated above, calibration parameters of 
GRBs and cosmological parameters are fitted simultane- 
ously, and measurement errors and systematic errors are 
both taken into account. We can see that the deviation 
from the cosmological constant at redshifts z > 0.5 turns 
out to be greater than the results in ||, such that the 
EOS of —1 lies almost at the edge of the 2a confidence 
interval. Though still consistent with the cosmological 
constant at 2a confidence level, such a deviation should 
be large enough to attract our attention. Of course, if 
this deviation is just an illustration of statistical errors 
due to the limitation of current observational data, it 
would be meaningless and should disappear with the in- 
crease of the observational data. While a comparison of 
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FIG. 1: Estimates of the uncorrelated dark energy EOS pa- 
rameters Wi. Top: uncorrelated dark energy parameters ver- 
sus redshift, in which the vertical errorbars correspond to la 
and 2a confidence levels of Wi and the horizontal errorbars 
span the corresponding redshift bins from which the contri- 
butions to Wi come most. Bottom: Window functions for 

Wi. 



the top plot of Figure Q] with Figure O for which GRBs 
are not included in constraining, shows that SNe la alone 
shift the dark energy EOS at redshifts z > 0.5 upwards 
from the cosmological constant and GRBs shift it a little 
more in the same direction. This means that both SNe la 
and GRBs prefer a dark energy EOS greater than —1 at 
redshifts z ^ 0.5. One can argue that the independence 
of SNe la and GRBs reduces the possibility of the devia- 
tion arising from statistical errors. Such a deviation from 
the cosmological constant, if confirmed, may be caused 
by the nature of the dark energy itself or some biasing 
systematic errors in the observational data that should 
be excluded. For the latter, we would need reconsider 
the process of calibrating SNe la and/or GRBs. It is no- 
table that the recent CfA3 addition of SN la samples has 
brought SN la cosmology to the point where systematic 
uncertainties dominate 2^, 23| . While the former is more 
exciting for possibly ruling out the cosmological constant 
as the dark energy. A close attention should be paid to 
this deviation with future observational data. 
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FIG. 2: Estimates of the uncorrelated dark energy EOS pa- 
rameters. Same as the top plot of Figure [1] except that GRBs 
are not included in constraining. 



To be more careful, we crosschecked our results by in- 
cluding GRBs in other different ways. For the results in 
Figure [U GRB data are not processed prior to being used 
to constrain the dark energy, i.e. the calibration of GRBs 
and constraining cosmological parameters are carried out 
simultaneously. In light of model-independent calibra- 
tions of GRBs in literatures 0, [13, II | , we performed the 
same analyses including instead the preprocessed GRB 



data from [ll[ and 24 1. In [ll[, GRB data are summa- 



rized by a set of model-independent distance measure- 
ments. These distance measurements can be used di- 
rectly to replace GRBs in constraining cosmological pa- 
rameters. In [iij ]. GRBs of redshift z < 1.4 are utilized 
to calibrate the luminosity relations based on a local re- 
gression estimate of distance moduli using the Union SN 
la sample 12], so the GRBs of redshifts z > 1.4, whose 
distance moduli are derived from the calibrated luminos- 
ity relations, can be used in the same way as SNe la. 
We present in Figure [3] the results of including GRBs in 
the above two ways. The bottom plot of Figure [21 for 
which the calibrated GRBs of redshifts z > 1.4 from [24| 
are used, is consistent with the results in Figure [TJ except 
that the constraints are slightly tighter such that the cos- 
mological constant has been ruled out at 2er confidence 
level at redshifts z > 0.5, which can be easily understood 
— some of the SNe la are used in both calibrating GRBs 
and constraining cosmological parameters. However, the 
top plot of Figure El for which GRBs are included by 
using the distance measurements from [TTJ] , is somewhat 
different from the results in Figure [U Comparing with 
the results without including GRBs (see Figure ^ , we 
can see that including these distance measurements does 
not change the result much. In fact, it was shown in [ill ] 
that these distance measurements shift best-fit parame- 
ter values towards the cosmological constant. Since the 
derivation of the distance measurements from GRBs in- 
volves quite a few intermediate steps and are carried out 
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FIG. 3: Estimates of the uncorrelated dark energy EOS pa- 
rameters. Same as the top plot of Figure [1] except that GRBs 
are included in model-independent ways. Top: GRBs are 
included by using the distance measurements from [ill ]. Bot- 
tom: GRBs are included by using calibrated GRBs of red- 
shifts z > 1.4 from [H|. 



through the Markov chain Monte Carlo method, it is ob- 
scure what has caused the difference. 

At the end, we would like to mention that in the above 
analyses we did not use the recent Constitution set of 
SNe la [23| and the BAO measurements presented in [25| . 
First, this is for the consistency of the data. Because, 
in 24 1, GRBs are calibrated with the Union set of SNe 
la 12j |. Second, there seem to be some tension in these 
data sets. The results derived using them are quite dif- 
ferent from the above. For the BAO measurements pre- 
sented in [2||, see the argument in [!§]• For SNe la, we 



noted that Union set prefers a Hubble parameter around 
70 km/s/Mpc. While Constitution set are derived by 
adding CfA3 SNe la to Union set using a Hubble pa- 
rameter of 65 km/s/Mpc. We wonder whether this will 
cause any problems of consistency or not. Anyway, we 
present in Figure [4] the results using Constitution set of 
SNe la [23j] and the BAO measurements from [iH, leav- 
ing the clarification of the differences between the data 
sets for the future. In spite of this, we can see from Fig- 
ure 0] that our conclusion on the dark energy EOS at 
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FIG. 4: Estimates of the uncorrelated dark energy EOS pa- 
rameters. Same as the top plot of Figure [1] except that Con- 
stitution set of SNe la from (2^] and the BAO measurements 
from ' 25(| are used instead. 



redshifts z > 0.5 are unaffected. See [26| for a discussion 
on the behaviour of the dark energy at low redshifts de- 
rived from Constitution set of SNe la [23J and the BAO 



measurements presented in [25 1. 

In summary, motivated by the fact that both SNe la 
and GRBs seem to prefer a dark energy EOS greater than 
— 1 at redshifts z > 0.5, we perform a careful investiga- 
tion on this situation, including more careful treatments 
of measurement errors of GRBs than previous studies on 
the evolution of the dark energy and crosscheck by using 
GRBs in different ways. We find that the deviation of 
dark energy from the cosmological constant at redshifts 
z > 0.5 is large enough that we should pay close atten- 
tion to it with future observational data. Such a devia- 
tion may arise from some biasing systematic errors in the 
handling of SNe la and/or GRBs or more interestingly 
from the nature of the dark energy itself. 
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In the published version of there are typos in the lu- 
minosity relations Eq. (16) and Eq. (17) and the typos 
are passed on to Eq. (21). However, the correct equations 
have been used in the calculations, so the results are un- 
affected. The typos have been corrected in third version 
on arXiv. 



